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ABSTRACT 

We analyze the impact of the Fermi non-detection of gamma-ray emission 
from clusters of galaxies on hadronic models for the origin of cluster radio ha- 
los. In hadronic models, the inelastic proton-proton collisions responsible for the 
production of the electron-positron population fueling the observed synchrotron 
radio emission yield a gamma-ray flux, from the decay of neutral pions, whose 
spectrum and normalization depend on the observed radio emissivity and on the 
cluster magnetic field. We thus infer lower limits on the average cluster mag- 
netic field in hadronic models from the Fermi gamma-ray limits. We also calcu- 
late the corresponding maximal energy density in cosmic rays and the minimal- 
guaranteed gamma-ray flux from hadronic radio-halo models. We find that the 
observationally most interesting cases correspond to clusters with large radio 
emissivities featuring soft spectra. Estimates of the central magnetic field values 
for those clusters are larger than, or close, to the largest magnetic field values in- 
ferred from Faraday rotation measures of clusters, placing tension on the hadronic 
origin of radio halos. In most cases, however, we find that the Fermi data do not 
per se rule out hadronic models for cluster radio halos as the expected gamma-ray 
flux can be pushed below the Fermi sensitivity for asymptotically large magnetic 
fields. We also find that cosmic rays do not contribute significantly to the cluster 
energy budget for nearby radio halo clusters. 
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INTRODUCTION 



The observation of diffuse, Mpc scale radio synchrotron emission in clusters of galaxies 
indicates the presence of both relativistic particles and large-scale magnetic fields in the 
intra-cluster medium (ICM). Morphologically, diffuse cluster radio emission is observed in 
two distinct varieties: radio relics and radio halos. Radio relics - polarized, asymmetric 
structures typically at the outskirts of clusters - may be associated to prima ry relativistic 



electron populations directl y accelerated at merger or accretion shocks (e.g. lEnsslin et al. 



19981 ; iKempner et al.l 12004 ). In contrast, radio halos are regular, un polarized, and dis 



tributed similarly to the th ermal X-ray emission from the ICM (see e.g. Ferrari et al.l 12008 
Feretti fc Giovanninill2008l for recent reviews). The large extent and regular appearance of 



radio halos is difficult to explain through a primary relativistic electron population as these 
particles are short lived, losing energy before they can diffuse through any significant distance 
in the cluster. Therefore, relativistic electrons should trace where they are produced. 

A couple of models have been proposed to explain the existence of radio halos. In the 
"hadronic" model, secondary cosmic-ray ele ctrons are generated in the collisions of rela - 
tivistic protons with nuclei in the ICM (e.g. iDennisonl Il980l ; iBlasi fc Colafrancescol Il999l ). 
Cosmic-ray protons, like primary electrons, can be accelerated in merger and accretion 
shocks, by AGN, or by supernovae; unlike primary electrons, though, cosmic-ray protons 
loose energy less efficiently, diffusing for cosmological time-scales, thus effectively building 
up and fairly homogeneously populating the cluster halo. A second possibility is that the 
relativistic electrons responsible for radio halo emission originate from lower e nergy electrons 
which are turbulently re-accelerated as a consequence of cluster mergers (e.g. 
200ll ; |Petrosian]l200lh . 



Brunetti et al. 



It has been argued that the h adronic model h as difficulty explaining the detailed prop 



erties of some radi o halos (see e.g. 



the radio emission ( jMarchegiani et al. 



Brunetti 



spectra observed for some halos (IBrunetti 



2007 



2004 , for a review) such as the radial extent of 



2009 



Macario et al.l 



the r adio spectru m with increasing freque ncy or radius (IFeretti et al.l 12004 ; iThierbach et al. 



Dqnnert et al.l 120101). th e steepness of the radio 



2010), and t he steepening of 



20031 ) (see however iKeshet &; Loebl (120101 ) ) . Other studies, however, have pointed out that 
hadronically generated cosmic-ray electrons emitting synchrotron radiation in the presence 
of a strong magnetic field (B > 3/xG) offer the most natural explanation of the strong cor- 
relation of radio halo power with cluster thermal X-ray luminosity and that merger induced 
enhancement and subsequent decay of the cluster magnetic field can explain th e associations 
of radio halos t o mer ging clusters (IKushnir et al.l 12009c IKeshet fc Loebl 120101 ) (see however 
Brunetti et al.l ( 120091 )). Models of cluster radio halos depend sensitively on the strength of 
the cluster magnetic field, on the cosmic-ray density, and on the distribution of both of these 
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quantities throughout the cluster volume, none of which is robustly measured at present. 
Some information on cluster magnetic fields can be gained through Faraday rotation mea- 
sures of radio galaxies in or behind clusters, though the translation of rotation measures 
to magnetic field strength depends sensitively on the unknown spa tial structure of cluster 
magnetic fields as we ll as on the gas density distribution (see e.g. iGovoni fc Ferettil 12004 ; 
Carilli fc Taylor! 120021 and references therein). 



In hadronic models of radio halos, the relativistic electrons responsible for the radio 
emission are generated through inelastic collisions of high-energy cosmic-ray protons with 
nuclei in the intra-cluster gas. These same collisions lead to gamma-ray emission through 
the production and subseque nt decays of neutral pions. Diffuse gamma-ray emission has yet 



to be detected from clus ters (IReimer et al.l 12003 



Aharonian et al. 2009; Aleksic et al. 2010 



Ackermann et al.l l2010al Jbl). placing constraints on the possible energy density in cosmic- 
ray protons and therefore on the density of relativistic electrons that may be produced 
from these. Recently, the Fermi-LAT Gamm a-ray Space Telescope ha s reported gamma- 
ray upper limits for a large sample of clusters (jAckermann et al.ll2010bi) which are roughl y 
an order of magnitude more stringent than previous EGRET limits (IReimer et al.ll2003l ). 
These limits imply that the energy density in cosmic-ray protons in nearby clusters is no 
more than a few percent of the thermal energy density, th ough the exact limits de pend on 
the assumed cosmic-ray spectral and spatial distributions (jAckermann et al.ll2010bf ). In the 
present analysis we qu a lify th e implications of the new gamma-ray observations for hadronic 



models (IReimer et al.l (120031 ) point out th at the EGRET data rule d out the gamma-ra y 
fluxes predicted in the hadronic models of lEnsslin et al . (1997) and Par &; Shavivl ( 119951 ) 



though these assume rather high cosmic ray densities. IReimer et al. rl2003h did not address 
the questions of the implications for the cluster magnetic fields or the minimal gamma-ray 
flux which we study here). Also, gamma-ray li mits from clusters of galaxies have significan t 
implications for particle dark matter searches ( jAckermann et al.ll2010al ; iDugger et al.ll2010l ). 



In this paper, we re-examine the hadronic scenario for the origin of radio halos in light of 
the tight upper limits placed by Fermi on the possible gamma-ray emission from clusters of 
galaxies. The gamma-ray non-detections effectively limit the possible density of secondary 
relativistic electrons and imply a minimum magnetic field strength if the observed radio 
halo fluxes are to be reproduced. In §2, we present our sample of radio halo clusters with 
Fermi gamma-ray upper limits. In §3, we use the observed radio fluxes and radio spectra 
under the assumption that radio halos are hadronically produced to predict the minimum 
possible gamma-ray emission and to extract the minimal magnetic field strength required 
not to exceed the Fermi gamma-ray limits. In §3.4, we compare the resulting magnetic field 
strengths required in the hadronic scenario to Faraday rotation estimates of cluster magnetic 
fields. Finally, we discuss our results and draw our conclusions in §4. 
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CLUSTER SAMPLE 



We consider the sample of nine, nearby radio halo clusters which have pub lished upper 

limits on their gamma-ray emission from the Fermi Gamma-ray Space Telescope (lAckermann et al 
2010bl ). These clusters represent the best cases to test the hadronic production of radio ha- 
los. We do not consider here cluster radio relics, since relics li kely have a differen t origin 
(e.g. primary electrons directly accelerated in merger shocks (e.g. lEnsslin et al.lll998l )). Most 
of the clusters in our sample host so-called giant radio halos with ~ Mpc scale radio emis- 
sion, although the halo in A2029 (and possibly in A2142) is classified as a mini-halo, i.e . 



a smaller diffuse radio sou rce found at the center of cool core clusters (IMurgia et al.l 12009 
Feretti k Giovanninill20081 ). 



The basic radio properties of our cluster sample (the radio flux at 1.4 GHz and the power 
law spectral index) collected from the literature are listed in Table 1, w hich also indicates 



the upp er limits on the gamma-ray flux in the 0.2-100 GeV band from lAckermann et al. 



(l2010bl ). In this work, we focus on whether the total radio flux and average radio spectrum 
of cluster radio halos can be reproduced in hadronic models given the Fermi gamma-ray 
upper limits. We will defer the consideration of the detailed spatial and spectral properties 
of radio halos, which may provide additional constraints, to future work. 



3. CONSTRAINTS ON HADRONIC MODELS 
3.1. Radio Halo Modeling 

We approximate the primary proton cosmic-ray population with a power-law 

dAL _ / E \~ ap 

with a low-energy cut-off set at E ~ 0.3 GeV, a value inferred from estimates of the minimal 
cosmic-ray energy provided by acceleration processes in clusters (e.g. shocks and supernovae) 
that does not insta ntly thermalizes due to Coulomb and hadronic energy losses (see e.g. 



Jubelgas et al.l 120081 ) 



Electron-positron (e^) production results from the decays of charged pions copiously 
created in inelastic scattering processes involving a primary cosmic-ray proton from the 
population in Eq. (CQ) and a target nucleus (essentially at rest) in the intra-cluster medium. 
The resulting secondary e ± population features an approximate injection spectrum of the 
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form (IMannheim &: Schlickeiserl Il994t iDermerl Il986l ) : 



dNe 
dE 



oc <7ppn N (f)N Q (f) 



( e y ap 
VGevy 



(2) 



where a pp ~ 32 mbarn indicates the inelastic proton-proton cross section at the energies of 
interest. 

The population of high-energy e^, which eventually will be responsible for the radio 
emission (as well as for inverse Compton emission potentially at energies as large as those 
detectable with Fermi), loses energy on short time scales corr esponding to a spatial diffusion 
length much smaller than the size of a typical cluster (see e.g. ISarazinlll999f). In other words. 



the effect of high-energy electron diffusion in clusters is always negligible (IBerezinsky et al. 



19971 ; IColafrancesco fc Blasilll998l ). and the steady-state solution for the equilibrium e ± pop- 
ulation results from accounting for energy losses alone, i.e. 



b(E e 



dE' 



dNe 
dE 



(3) 



wher e b(E e ) indicates the energy loss term in the diffusion-loss equation (see e.g. IColafrancesco et al 
20061 ). We include in the function b(E e ) all relevant energy loss processes, including in- 
verse Compton scattering of CMB photons, synchrotron emission (assuming an average 
constant magnetic field in the region of interest), bremsstrahlung (u sing as an input the 
appropriate beta models for the gas density, from I Chen et al.l (120071 )) and Coulomb scat- 
tering. After having established the steady-state equilibrium e ± solution, we calculate the 



sy nchrotron, inverse Compton and bremsstrahlung emi ssion following the 



m 



Colafrancesco et al.l (120061 ) and recently reviewed in lProfumo fc Ulliol (120101 ) 



ormalism outlined 



The dominant gamma-ray emission process in the hadronic scenario is via the two- 
photon decay of neutral pions, which are produced with approximately the same multiplicity 
as their charged counterparts. To model the resu lt ing g amma-ray spe ctrum we follow the 
analyt ic approximation to the results of IDermerl (119861 ) proposed in iPfrommer fc EnBlin 
(120041 ) . We verified that the relative normalization in the injection spectrum of gamma-rays 
from neutral pion decays vers us electrons from charg ed pion decays follows the approximate 
relation, at E ^> m^, (see e.g. iDolag fc EnBlinll2000l ) 



dl% 
~dE 



(E) 



d^ 



16a„ - 8 dE 



(4) 



-6 - 



3.2. Limits on Cluster Magnetic Fields, Gamma-ray Emission, and 

Cosmic-Ray Energy Densities 

Fixing a value for the average magnetic field across a given radio halo, the observed 
radio flux yields the normalization to the primary cosmic-ray proton density In turn, this 
produces a unique prediction (with the present assumptions on the morphology of cosmic-ray 
sources, and assuming a constant magnetic field intensity across the region of interest - see 
the next section for the effect of a non-trivial magnetic field morphology) for the gamma-ray 
flux from neutral pion decay, as well as for the sub-dominant bremsstrahlung and inverse 
Compton emission processes from the secondary electron-positron population. Comparing 
this prediction to the Fermi limits on the gamma-ray flux from clusters exhibiting a radio 
halo, we can infer the minimal value of the cluster average magnetic field compatible with 
gamma-ray observations: since the radio emissivity is proportional to the square of the 
magnetic field intensity for a given electron density distribution, the smaller the magnetic 
field the larger the relative cosmic-ray proton primary population normalization, and hence 
the larger the gamma-ray flux from neutral pion decay. 

We list in Table 2 our results for the constraints on hadronic models for cluster radio 
halos from Fermi observations. A crucial input parameter is the slope of the spectrum of 
the primary cosmic-ray protons, i.e. the parameter a p in Eq. (TTJ. We fix this parameter to 
match the slope of the radio emission at a given frequency, as listed in the fourth column 
of Table 1. For a power-law steady-state electron distributio n with a spectral index ay , the 



radio emissivity has a corresponding spectral index (see e.g. iRybicki &: Lightmanlll986l ) 

«H = («e - l)/2. 

For the present case, in the limit of high-energy, b(E) ~ E 2 and for the injection spectrum 
given in Eq. (T5]) one thus has 

n ^TT^ [ (ET ap dE' ~ 
°{E) J E 

which therefore implies 

~ Qtp/2. 

Of course, the radio spectrum depends on the assumed magnetic field as well as, potentially, 
on other energy loss terms contributing in b(E) which are not quadratic in energy, such as 
bremsstrahlung, so the above is an approximation. Nonetheless, this is a good guideline to 
guess the input a p for which one can obtain the correct radio spectrum, and we then iterate 
considering the full energy loss equation to find a p . Comparing the values of cir from Table 
1 and of a p from Table 2, we do find a good confirmation of the theoretical approximation. 
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Ackermann et al.l (l2010bf ) report Fermi gamma-ray upper limits for a large sample of 
clusters in four energy bands in their Table 2. Having fixed the appropriate cosmic-ray spec- 
trum, we calculate the minimal average magnetic field such that the predicted gamma-ray 
flux does not exceed the limits quoted in any one of the four Fermi bands. The resulting 
minimum average magnetic fields and the Fermi bands giving the most stringent constraints 
for each individual cluster are listed in the third and fourth columns of our Table 2, respec- 
tively. Notice that the values of the magnetic fields we quote here are "average" in the sense 
that we assume a single magnetic field value for the radio halo. 

If cluster radio halos originate from hadronic cosmic-ray processes, the larger the average 
cluster magnetic field, the smaller the predicted gamma-ray flux. The latter, however, is 
bounded from below: we provide here a calculation of this minimal gamma-ray flux expected 
in the context of hadronic radio halo models. In the limit in which the magnetic field 
intensity is much larger than the equivalent magnetic field energy density in the CMB, 
B 2 ^> -E>cmb> where -Bcmb = 3.2(1 + z) 2 fiG, synchrotron dominates the energy losses of high- 
energy electrons, and n e ~ B~ 2 . Since the radio emissivity is proportional to the magnetic 
field intensity squared and directly proportional to n e , in the limit B 2 ^> -B^mb the radio 
emissivity is essentially independent of the magnetic field intensity. Correspondingly, in this 
limit one gets the smallest possible value for the primary cosmic-ray proton density, and 
hence the minimal 7-ray flux in the context of the hadronic scenario for cluster radio halos. 
The fifth column of Table 2 indicates how far the current Fermi limits are from this minimal 
( "guaranteed 1 ) flux of 7 rays, for the best Fermi energy band (indicated in the column to the 
left). We find that the current Fermi sensitivity is anywhere between 0.1% and almost 10% 
of this guaranteed 7-ray flux in the hadronic scenario, where the least promising clusters 
have the hardest cosmic-ray proton spectra (a p < 3) while the most promising clusters have 
the softest spectra (a p > 3). This is easily understood, since soft cosmic-ray spectra produce 
copious soft neutral pions that dump 7-rays within the low-energy range accessible to Fermi, 
while they hardly produce any electrons-positrons energetic enough to radiate synchrotron 
radio emission at the required level. This implies that to achieve the measured radio flux 
the normalization of the cosmic-ray primaries must be tuned to large values, yielding, in 
turn, very large 7-ray fluxes. We also indicate, in the sixth column of Table 2, the integrated 
minimum 7-ray flux for the hadronic model across the entire Fermi energy range (0.2-100 
GeV), to be compared with the sixth column of Table 1 (showing, in the same energy range 
and with the same units, the current limits). 

The last column of Table 2 shows our results for the maximal value of the ratio €crI e th of 
the cosmic-ray proton (kinetic) energy density 6cr versus the cluster thermal energy density 
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compatible with the Fermi limits. Here, n ejt h indicates the number density of thermalized 
electrons, 

1 - lx T 



cL, 



1 + 



-He 



i-lx 



He 



is the number of particles per electron in the ICM, and we set the 4 He mass fraction to its 
primordial value Xn e — 0.24. Finally, we assume that the ICM temperature T is constant 
with radius. To calculate the cosmic-ray energy density, we use the cosmic-ray primary 
spectrum outlined above (featuring a cut-off at 0.3 GeV for all clusters), with a p set to the 
values indicated in the second column of Table 2. We assume that the spatial distribution of 
the primary cosmic-ray protons follows the intra-cluster thermal gas density (i.e.: iVo(r) oc 
n e,th{^)), which we model using an isot hermal beta - mode l. We employ the temperature values 
and beta-model parameters listed in I Chen et al.l (120071 ). Notice that with our assumption 
of cosmic-ray proton sources following the ICM gas density, and of negligible cosmic-ray 
diffusion, the ratio ecnl^th is a constant in space. We also note that the figures we quote 
in the last column of Table 2 depend on our assumption that the cosmic ray proton density 
scales like the thermal gas density. Relaxing this assumption and allowing for a relative 
bias in the two density distributions might produce a scatt er by a factor of a fe w in the 
constraints on the relative cosmic ray energy density, see e.g. iJeltema et al.l ( 120091 ) . 



We find upper limits on the ratio ecn/tth ranging between 0.01 to over 1.0 for the radio 
halo clusters in our sample. For the nearest radio hal o clusters, cosmic rays c ontribute no 
more than 1-5% of the thermal energy density (see also lAckermann et al.ll2010bl ). The values 
listed in column 7 of Table 2 are "maximal" in that they refer to the minimal value of the 
cluster magnetic field inferred from the radio emission: increasing the magnetic field sup- 
presses the cosmic-ray normalization while not affecting substantially the associated spectral 
shape. This is illustrated visually, for the case of three clusters (Coma, A1914 and A2163) 
in Fig. 1, where we plot ecnl^-th as a function of the average cluster magnetic field. We indi- 
cate the impact of the Fermi limits, which constrain the magnetic fields to be above certain 
values, and therefore the energy density ratio to be below certain corresponding levels. The 
figure illustrates in the hadronic model clusters with a soft cosmic-ray spectrum (A1914) 
feature very large relative cosmic-ray energy densities, which, for reasonable magnetic field - 



values, are typically more than 10% of the thermal energy density (see also iBrunettil 12009 



1. 



Macario et al .1 120 10 ). The opposite is true for clusters with hard cosmic-ray spectra. No- 
tice that the relative normalization of the curves depends sensitively on the cluster radio 
emissivity, temperature and gas density profile. 

It shoul d be noted that ra dio halos are often observed to be associated to merging 
clusters (e.g. lCassano et aU2010[ ). and several of the clusters in our sample exhibit non-trivial 
spatial variations in the ICM density and temperature. In the hadronic model, gamma- 
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ray emission and relativistic electrons result from the same proton-proton collisions, so the 
calculation of the minimum gamma-ray emission does not depend assumptions of the ICM 
distribution or temperature nor does the derivation of the minimum magnetic field. We have 
made the simplifying assumption of a constant average magnetic field (note the minimum 
gamma-ray emission does not depend on the magnetic field), which is likely not the case, 
but without a resolved gamma-ray or non-thermal X-ray detection, spatial variation in the 
magnetic field and cosmic ray electron population cannot be separated; as clusters are more 
or less point sources for Fermi, the gamma-ray limits give us only a single average constraint 
on the cosmic ray population. Below we consider a couple of simple magnetic field radial 
profiles. Variations in the ICM density and temperature will effect the estimate of the cosmic 
ray energy density and thermal energy density (though here are average temperature should 
be reasonably accurate). Under the assumption made above that the cosmic ray density 
traces the thermal gas density, the effect is minimal, but as already mentioned biases in the 
location of the cosmic ray population relative to the thermal gas will effect the constraints 
(at the level of a factor of a few for cosmi c ray distributions motivated by simulations (e.g. 



Pfrommer et all 120081 : 1 Jeltema et al.ll2009j )). 



3.3. Estimates for the Central Magnetic Field 

The results we quote above refer to the simplified scenario where the cluster magnetic 
field is constant throughout the cluster. The limits on the magnetic field quoted in Table 2 
obviously depend on this assumption. In this section, we consider a more realistic radially 
varying magnetic field and, in particular, we extract estimates for the central magnetic field 
value. In the next section, we then compare our minimum required magnetic fields in the 
hadronic model to Faraday rotation measurement determinations of cluster magnetic fields 
and ask whether the needed magnetic fields actually exist in clusters. 

Both MHD simulations and Faraday rotation measures, for the few clusters with rotation 
measures at multiple distances, i mply that cluster magnetic fields scale with the ICM gas 



2009 



densi t y and decrease with r adius ( iDolag et al.ll200ll ; iGovoni fc Ferettil 12004 iDonnert et al. 



Bonafede et al.ll2010l ). Denoting with F(r) the spatial distribution of the thermal 



cluster gas as given by a beta-model with core radius r c , i.e. 



n N {r) 
n N (0) 



F(r) 




and assuming that the magnetic field scales with the gas density, we employ a phenomeno- 
logical radial dependence of the magnetic field strength - to be motivated below - of the 



- 10 - 



type: 

B(r) = B c [F(r)]". 

In this case, the central value of the magnetic field B c can be expressed as a function of the 
average magnetic field Bq (which we constrain in the previous section) as: 

R2 = 2 / r "r 2 [F(r)] 2 dr 
c °/ rh r 2 [F(r)] 2 + 2 Mr' {) 

where indicates the radial size of the radio halo. It is straightforward to generalize the 
equation above to any magnetic field spatial distribution. 

We quote the constraints on the central magnetic field for two cases, rj = 0.5 and rj — 1, 
for the radio halo clusters under consideration here in the last two columns of Table 3. These 
values for rj effectively bound the values found in simulations and Faraday rotation measure 



studies. For example, rj = 0.5 provides a good match to the radial dependence of 



rotation measures observed for the Coma cluster in the recent study of Bonaf ede et al. 



' arada y 

fcoioh . 



while simulations find a value of 77 close to one (IDolag et al.ll2001l ; iDonnert et al.ll20090. and 
77 = .9 is a good match to the rotation measures observed for the cluster A119 (IDolag et al. 



20011 ). A magnetic field frozen in to the plasma would give instead an intermediate value, 
77 = 2/3. For reference, we note in the second and third columns of Table 3 the radio halo 
size we consider and the minimal average magnetic field, -Bo, determined from the Fermi 
data in the analysis in the previous section. 



3.4. Comparison to Faraday Rotation Measures 



The only cluster in our sample whose magnetic field has been s tudied in detail 



Faraday rotation measures is A1656 , also known as the Co ma cluster (IBonafede et al.ll2010 



using 



Kim et al. 1990 : Feretti et al. 1995 ). The recent study of Bonafede et al. fcoid ) analyzed 
polarization data for seven radio sources within the Coma cluster at varying radii to derive the 
radial profile of the magnetic field strength. They find a best fit central magnetic field value 
of B c = 4.7/iG for a radial dependence compared to the thermal gas of rj = 0.5, consistent 
with but close to the value of 3.9/iG for t he same 77 that we infe r for the minimum central 
field in the hadronic scenario. Similarly, IBonafede et al.l ( 120101 ) find an average magnetic 
field over the cluster volume of ~ 2/iG close to our minimum derived value of B = 1.7/iG, 
and for rj = 1 our central magnetic field of 7.4/zG is only marginally consistent with the 
rotation measure data (~ 1% confidence level). Therefore, while the hadronic generation of 
the radio halo in the Coma cluster is still consistent with Faraday rotation measure estimates 
of its magnetic field, any improvement in the gamma-ray upper limits will make the needed 
magnetic field in hadronic models too high compared to Faraday rotation measure estimates. 
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Among our cluster sample, A1914 which hosts a very steep spectrum radio halo, notably 
exhibits a very high required magnetic field if the radio halo is of hadronic origin. For 
this cluster, we derive a minimum average magnetic field of 7.4/zG and a minimum central 
magnetic field in the range 21 — 40/iG. While Faraday rotation measure estimates of the 
magnetic field are not available for this cluster, these values are much higher than the 
typical few /iG magnetic fields found in rotation measure studies of non-cooling flow clusters 
like A1914 (and of clusters hosting giant radio halos in general which tend to be mergers); 
Faraday rotation studies only imply such high values of the magnetic field for observations 
of radio galaxies in the centers of c ool core clusters (in some cases hosting radio mini-halos ) 



and not over such large scales (e.g. Ferrari et al.ll2008l ; ICarilli fc Taylorl l2002r IClarkd 120041 ) . 



The requirement of such a large average magnetic field in A1914, therefore, sheds doubts on 
a hadronic origin for this radio halo. In addition to Coma and A1914, the two clusters A2256 
and A2319 also have minimum required magnetic fields of a few /zG which, while consistent 
with the typical magnetic fields derived from Faraday rotation measures, will similarly soon 
be in tension with these if clusters continue to be undetected with deeper Fermi gamma-ray 
observations, or perhaps with new Faraday rotation measure data. A2319, for example, has 



an ea rly Faraday rotation measure estimate of its magnetic field strength of 2/iG ( jVallee et al. 



1983). 



We warn the reader that a few assumptions (including the functional dependence of 
the magnetic field power spectrum and the magnetic field radial distribution) come into 
play w hen translating Fara day rotation measures into de-projected magnetic field values, 
see e.g. iMurgia et al.l (120041 ) and lBonafede et al.l (120101 ) . Finally, we note that equipartition 
estimates of the average magnetic fields in radio halo clusters are typically much smaller than 
those derived from Faraday rotation measures, being in the range of 0.1 — 1/iG ( iFerrari et al. 



20081 ) , and are smaller than the minimum values derived here for many of the clusters in our 
sample. 



4. DISCUSSION 

In this paper, we analyze how early Fermi constraints on the gamma-ray flux from clus- 
ters of galaxies impact hadronic models for clusters radio halos. We consider all clusters 
with both a radio-halo detection and for which Fermi gamma-ray upper limits are available. 
We tune the cosmic-ray proton spectrum responsible, in hadronic models, for the subse- 
quent production of the non-thermal electron-positron population that fuels the synchrotron 
emission observed at radio frequencies, to match the slope and intensity of the detected 
radio emission. We then calculate the smallest possible value for the average magnetic field 
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in the radio halo region compatible with both the total observed radio halo flux and with 
the Fermi limits on the gamma-ray emission which would inevitably result from the same 
inelastic proton-proton collisions producing the electron-positron population. We note that 
the resulting magnetic field constraints are especially interesting for those radio halos with 
a soft radio spectra. For clusters such as A1914 or A2256, we find that in order to avoid 
over-producing gamma rays at Fermi energies, average magnetic fields at least on the order 
of several micro-Gauss are needed. Using two phenomenological models for the magnetic 
field radial dependence, the required central magnetic fields are a few /xG for several clus- 
ters and over 20 /iG for A1914. We find that in a few cases the implied magnetic fields in 
hadronic models are in excess of, or close to, the largest cluster magnetic field values obtained 
with Faraday rotation measures, placing tension on the hadronic origin of radio halos. The 
large values of the average magnetic fields we infer also imply a large non-thermal pressure, 
well beyond equipartition, that would have dramatic implications for the assumption of hy- 
drostatic equilibrium and for cluster mass estimates and their use for cosmological studies. 
Further information may be gained by considering spatial variation in the observed radio 
flux and spectrum for specific selected clusters with the required wealth of radio data, and 
we defer this question to a forthcoming study (Jeltema et al, in prep.). 

We also estimate the ratio of the energy density in cosmic rays over the cluster thermal 
energy density - a quantity that is maximal for the minimal value of the magnetic field we 
infer. We find that this ratio varies between a percent to values of order unity, with larger 
values associated to both very hard (cosmic-ray spectral index a p ~ 2) or very soft (a p > 3) 
cosmic-ray - and hence radio - spectra. For at least a few nearby radio halo clusters it 
is cl ear that cosmic rays do not contribute significantly to the cluster energy budget (see 



also lAckermann et al.l l2010bl ). The minimal gamma-ray emission expected in the context 
of hadronic radio-halo models, corresponding to asymptotically large values of the cluster 
average magnetic field, is also derived. This minimum emission is typically between one to 
three orders of magnitude below the current Fermi sensitivity. The most promising cases are 
those corresponding to clusters with soft radio spectra, for which the needed improvement 
in sensitivity to detect a gamma-ray emission associated to a hadronic radio halo for any 
magnetic field is on the order of 10. 

In summary, we showed that while hadronic models for radio halos are not ruled out 
by Fermi gamma-ray limits on galaxy clusters, the required cluster magnetic fields are in 
some cases high and possibly in tension with other observations. This result reenforces how 
the high-energy end of the electromagnetic spectrum provides an important and crucial tool 
to understand the nature of non-thermal phenomena in the largest bound structures in the 
universe. 
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Table 1. Cluster Data 



Cluster 


z 


1.4 GHz Flux 




Ref. 


7-ray Upper Limit 






Density (Jy) 






(10~ 9 ph cm- 2 s" 1 , 0.2-100 GeV) 


A1656 


0.023 


0.64 


1.35 


1 


4.58 


A1914 


0.171 


0.064 


1.88 


2, 3 


1.19 


A2029 


0.076 


0.022 


1.35 


4,5 


3.28 


A2142 


0.090 


0.0183 


1.5 


6 


2.82 


A2163 


0.203 


0.155 


1.18 


7,8 


5.51 


A2256 


0.060 


0.103 


1.61 


9 


1.96 


A2319 


0.056 


0.212 


1.28 


4, 10 


0.75 


A2744 


0.308 


0.0571 


1.0 


11, 12 


2.49 


IE 0657-56 


0.296 


0.078 


1.3 


13 


2.75 



Note, 
references: 



Ga mma-ray upper l i mits a re ta ken from lAckermann et a" 



'D iThierbach et all (120031). (2) Bacchi et all (120031 /31 



(2010b). Radio 



Giovannini et al. 



(hoogh . (4) iMurria et all (boOfl h . (5) Isiee fc Siegmanl (Il983 l)~( 6) iGiovannini fc Feretti 



(I2000f). (7) 



Feretti et all (120041) . (8) 



Feretti et a 



Kempner fc SarazirJ (boOlh . (11) bovoni et all fcoOll ). (12) lorru et all (120071 ) . (13) 



jgOOlh . (9 ) ferentienj d2008l) . (10) 



Liang et all (120001 ) . 
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Table 2. Hadronic Model Results 



Cluster 


aip 


-Bo, mm 


Fermi Band 


Ratio of F 7jTnm 


Minimal 7-ray Flux 


Maximum 






( M G) 


(GeV) 


to Fermi up. lim. 


(10~ 9 ph cm -2 s _1 ) 


€ Cr/ e th 


A1656 


2.70 


1.70 


1-10 


0.040 


0.186 


0.018 


A1914 


3.76 


7.40 


0.2-1 


0.077 


0.181 


0.684 


A2029 


2.70 


0.37 


0.2-1 


0.001 


0.006 


0.052 


A2142 


3.01 


0.65 


0.2-1 


0.002 


0.010 


0.091 


A2163 


2.05 


0.37 


0.2-1 


0.004 


0.011 


0.131 


A2256 


3.25 


2.70 


0.2-1 


0.037 


0.099 


0.138 


A2319 


2.57 


1.40 


0.2-1 


0.040 


0.046 


0.010 


A2744 


2.01 


0.26 


1-10 


0.004 


0.004 


1.396 


IE 0657-56 


2.60 


0.77 


1-10 


0.005 


0.018 


0.508 



Note. - Column 2 lists the cosmic-ray proton injection spectral index implied by the radio 
spectrum for a magnetic field strength of B = 1/xG. Column 3 indicates the minimum average 
magnetic field required to produce the observed radio halo flux with secondary electrons without 
overproducing the gamma-ray flux compared to the Fermi limits. Column 5 gives the ratio of the 
minimum gamma-ray flux for any magnetic field value, corresponding to very large magnetic field 
values, in the hadronic model compared to the Fermi limits. Col u mns 3 and 5 are calculated for 
the most constraining Fermi energy band from lAckermann et all (l2010bl ). which is in turn noted 
in Column 4. Column 6, instead, gives the minimal integrated gamma-ray flux between 0.2 and 
100 GeV predicted in the hadronic scenario. The final column lists the maximum allowed ratio of 
the energy density in cosmic rays to the cluster thermal energy density implied by the Fermi limits 
given the proton spectral index implied by the radio spectrum and assuming that the cosmic-ray 
density traces the thermal gas density. 
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Fig. 1. — Ratio of the energy density in cosmic rays to the thermal energy density as a 
function of the average magnetic field, -Bo, for Coma (black solid line), A1914 (red dashed 
line) and A2163 (dot-dashed blue line). Also indicated on each curve are the minimum 
average magnetic field and corresponding maximum allowed energy density ratio needed to 
produce the radio halo flux without overproducing the Fermi gamma-ray upper limits. 
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Table 3. Magnetic Field Model Results 



Cluster 




Bmin 


B c , X] = 1.0 


B c , rj = 0.5 




(Mpc) 


(jmG) 


OiG) 


(jjG) 


A1656 


0.83 


1.70 


7.24 


3.91 


A1914 


1.04 


7.40 


40.1 


21.1 


A2029 


0.25 


0.37 


1.80 


0.93 


A2142 


0.19 


0.65 


1.51 


1.01 


A2163 


2.28 


0.37 


1.92 


1.03 


A2256 


0.81 


2.70 


8.48 


5.12 


A2319 


1.02 


1.40 


7.62 


3.82 


A2744 


1.89 


0.26 


1.15 


0.65 


IE 0657-56 


2.10 


0.77 


5.73 


2.78 



Note. 



Column 2 indicates the size of the radio halo 



we employ (IGiovannini et al.l l2009t iMurgia et al.l l2009t 
Giovannini fc Ferettillioooh . Column 3 lists the minimum 



average magnetic field required to produce the observed 
radio halo flux with secondary electrons without over- 
producing the gamma-ray flux compared to the Fermi 
limits. Column 4 and 5 list the predictions for the mini- 
mal value of the cluster central magnetic field, for 77 = 1.0 
and 7] = 0.5, where 77 indicates the bias in the magnetic 
field radial distribution with respect to the the cluster 
gas density distribution modeled beta-model. 



